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Abstract 

We describe the various techniques developed in the Fermilab Wideband Experi- 
ments, E687 and FOCUS, to reconstruct long-hved states. The techniques all involve 
modifications to standard tracking techniques and are useful to report for future ex- 
periments. 
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1 Introduction 



FOCUS is a photoproduction experiment located at Fermilab which has been 
configured to investigate the production and decay of charmed particles. As 
charmed particles typically decay to states containing strange particles, tech- 
niques have been developed to reconstruct long-lived mesons and baryons con- 
taining strange quarks. The techniques employ a full range of information from 
a multipurpose spectrometer including partially reconstructed tracks and kine- 
matic constraints, and utilize excellent charged particle identification using 
Cerenkov detectors, efficient muon identification, and the ability to recon- 
struct neutral hadronic energy. 

In this paper we will briefiy describe the spectrometer in Section 2. In Sec- 
tion 3 we will discuss how the tracking is performed in the silicon microstrip 
detector and in the multiwire proportional chambers and how the individual 
segments are linked to form tracks. In Section 4 we describe our neutral Vee 
reconstruction techniques. In Section 5 we present our algorithms for identi- 
fying E decays where one decay product goes undetected (Kinks). In Section 
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6 we present our methods for combining of the Vee and Kink algorithms to 
find S~ and fl~ decays. 



2 Spectrometer Components and Layout 



The spectrometer shown in Fig. 1 consists of two dipole magnets, Ml and 
M2; a 12 plane silicon microstrip array; five stations of multiwire propor- 
tional chambers; 300 threshold Cerenkov cells arranged in three gas boxes; 
two electromagnetic calorimeters for photon, electron, and 7r° reconstruction; 
two muon detectors; and a hadronic calorimeter for triggering and neutral 
hadron reconstruction. There are 10 meters along the spectrometer in which 
Vees can can be reconstructed. There is also a large charged particle accep- 
tance covering the entire forward hemisphere. The FOCUS spectrometer is an 
upgraded version of the previous experiment E687. A detailed description of 
the E687 spectrometer and its performance may be found in Reference 1. 
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Fig. 1. A schematic drawing of the FOCUS spectrometer. The inset displays the 
segmented targets, the embedded target silicon, the trigger counters, and the 12 
plane silicon tracking array. The spectrometer is about 30 meters long. 

The two large aperture magnets are used to momentum analyze charged 
tracks. The magnets are oriented to bend particles in the vertical (y-axis) 
direction with Ml and M2 bending in opposite directions. Positive particles 
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are bent upwards in Ml. The particular arrangement of magnet positions and 
orientation was chosen for its unique effect on event topology. There is a large 
background of e^e~ pairs coming from beam photon conversions in the ex- 
perimental target. Pairs are produced with little transverse momentum with 
respect to the beam direction {z direction) and with a transverse profile com- 
parable to the beam size (about 1 cm in x and y). The first magnet bends the 
electrons and positrons in y, creating a vertical swath. The lowest energy pairs 
strike the inside of Ml or the upstream face of M2, while the remainder pass 
through the M2 aperture and are bent back towards the beam axis. The beam 
profile is reconstituted at the end of the spectrometer, with some smearing due 
to energy loss via bremsstrahlung from material throughout the spectrometer. 

The target is composed of segmented BeO slabs with silicon microstrip dou- 
blets after the first two targets and after the last two targets. For historical 
reasons these four silicon microstrip planes are not used directly in the track 

reconstruction, however they are used to extend the microstrip tracks to the 
production and secondary vertices. By making a measurement close to the 
vertex we are able to improve both the position and the angular resolution of 
the track. 



3 Standeird Tracking Devices and Algorithms 
3.1 Microstrip Detector 

The high resolution tracking of charged particles provided by the Silicon mi- 
croStrip Detector, or the SSD, is essential to the reconstruction of charmed 
particle decay vertices and in separating these vertices from the production 
vertex. 

The SSD consists of four physically separated stations with three planes per 
station. The three planes of each station are oriented to i, j, and k coordinates 
of -135, -45, and -90 degrees with respect to the horizontal {x axis). The mi- 
crostrip detector and target layout is shown in Fig. 2 and a detailed summary 
of the active area, high resolution region, strip pitch, and number of channels 
for each station is given in Table 1. 

The innermost section of each plane, covering the region where tracks pass 
most closely to each other, has twice the resolution of the outer section. In 

addition, the most upstream station has a resolution twice as good as the other 
stations. Pulse heights are read out for all 8,256 strips. The overall detection 
efficiency of each plane is better than 99% even including the non-functioning 
strips and broken electronics channels. A more complete description of this 
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detector and its performance can be found in Reference 2. 
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Fig. 2. A layout drawing of tlie silicon microstrip detector together with the BeO 
production targets and embedded target silicon strip planes. 



The microstrip tracking algorithm is based on projection finding using hits in 
three separate views. A hit is a cluster of one to three adjacent strips depending 
on the summed pulse heights in the strips (less than a 1.5 minimum ionizing 
pulse height) and on the number of neighboring strips with pulse heights above 
threshold. As momentum information is not yet available multiple Coulomb 
scattering effects are not considered. Projections are found using very loose 
selection criteria. At least three hits per view (out of a possible four) are 
required and sharing of hits is permitted. Hits in the last three stations that 
are already assigned to a projection having hits in all four stations cannot 
be reused for a new projection containing only three hits. Projections are 
combined to form tracks if they match in space with a global pcr degree 
of freedom less than 8. Tracks sharing projections are arbitrated according to 
their x^- The process is performed in a fully symmetric way with respect to 
the views. 



The pulse height information on each hit is used in three ways. First, the 
information allows for the separation of two overlapping hits in the pattern 
recognition. Second, the pulse height information is combined from all hits in 
the SSD in the track to determine whether the track is consistent with the 
passage of a singly charged track through the planes. By requiring more than 
1.5 times a standard singly charged pulse height we can identify e+e~ pairs. 
Third, and most importantly, we use the charge sharing between neighboring 
strips to obtain an improved position resolution. 
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Table 1 

Properties of the Silicon Microstrip Detector. The first three lines are the internal 
relative z locations of the silicon microstrip planes. 



Property 


I station 


II station 


III station 


IV station 


1** detector 


-0.5 cm 


5.5 cm 


11.5 cm 


23.5 cm 


2"*^ detector 


0.0 cm 


6.0 cm 


12.0 cm 


24.0 cm 


3^'^ detector 


0.5 cm 


6.5 cm 


12.5 cm 


24.5 cm 


active area 


2.5 X 3.5 cm^ 


5x5 cm^ 


5x5 cm^ 


5x5 cm^ 


high res area 


1.0 X 3.5 cm^ 


2x5 cm^ 


2x5 cm^ 


2x5 cm^ 


strip pitch 


25 /xm, 50 /xm 


50 fim, 100 /Ltm 


50 fim, 100 /Ltm 


50 fj,m, 100 fxm. 


# of channels 


688 X 3 


688 x 3 


688 x 3 


688 x 3 



3.2 Multiwire Proportional Chambers (MWPCs) 



Five stations of Multiwire Proportional Chambers are used to track charged 

particles in the main spectrometer. The first three chamber stations, labelled 
PO, PI, and P2, are located between the two analyzing magnets. Ml and M2. 
The other two chamber stations, P3 and P4, are located downstream of the 
second magnet, M2. This arrangement allows for two momentum measure- 
ments for higher energy tracks which pass through both magnets as well as 
providing momentum information for wide angle tracks which are not accepted 
by the M2 aperture. 

Each of the chamber stations consist of four planes of wires, measuring y, u, v, 
and X positions. The X-view wires, running vertically, measure the horizontal 
(non-bend view) position. The U-view and V-view wires are oriented at ± 11.3 
degrees with respect to the Y-view, an arrangement which is used to resolve 
ambiguities and to provide better momentum resolution. 

The tracking algorithm for the multiwire proportional chambers is used to 
reconstruct spectrometer tracks which have the first hit in PO (the most up- 
stream chamber) and which extend into at least one additional chamber. As 
with microstrip tracks, a projection finding algorithm is used. The X-view is 
special because it is the only MWPC view in which the charged particles are 
not bent. 

First a microstrip track is extended into the non-bend X-view of the chambers. 
A search is made for MWPC hits which match this "seed" projection. Next, 
projections are formed from hits in the U-view, V-view, and Y-view. The 
projections in all four views are combined to form tracks. Further X-view 
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projections are formed using unused MWPC hits where no seed track exists 
and these projections are also matched with unused U-vicw, V-vicw. and Y- 
view projections to reconstruct additional tracks. Tracks are required to extend 
through at least the first three MWPCs. A least squares fit is performed on 
all candidate tracks. The fit parameters are the intercepts and slopes of the 
track in the xz and yz planes. If a track passes through just the first three 
MWPCs, then it is referred to as a stub. For tracks passing through all five 
chamber stations (called 5-chamber tracks or just tracks), the change in the y 
slope between the track segments upstream and downstream of M2, i.e. the 
bend angle, is an additional parameter. 

An approximate momentum is assigned to 5-chamber tracks by using the bend 
angle in M2 and the sudden bend approximation. Various magnetic corrections 
are applied to this category of tracks. Fringe fields in Ml and M2, extending 
beyond the magnetic poles, and off-field components of the magnetic field 
(By and Bz) are treated by higher order corrections to the linear fit. The 
corrections are applied by iterating each 5-chamber track through a complete 
fit which inchidcs the magnetic corrections. The track momentum is adjusted 
appropriately after each iteration. 



3.3 Linking 



Once tracks are reconstructed in the silicon microstrip and MWPC systems, 
they must be linked together to determine which microstrip track is associated 
to which MWPC track. This is accomplished by extrapolating each microstrip 
track to the center of Ml and searching for MWPC tracks which, when ex- 
trapolated back to this same point in z, "match" with it. This matching is 
performed by comparing the MWPC track's extrapolated x slope and x inter- 
cept at the center of Ml with the same quantities for the microstrip track in 
question. Because the magnetic field is in the x direction there is no (or little) 
bending of charged tracks in this view. Thus, for a MWPC track and a mi- 
crostrip track created by the same particle these two x quantities should agree. 
(This match is accomplished after a correction is made for weak focusing in 
the non-bend view.) A global least squares fit using all the hit information 
from the microstrip segment and the MWPC segment is then performed on 
each candidate. Multiple MWPC segments linked to the same microstrip seg- 
ment are arbitrated on the basis of the per degree of freedom of the fit. 
A maximum of two MWPC tracks are allowed to be linked to the same mi- 
crostrip track. This decision is made because e~^e~ pairs from beam photon 
conversions frequently reconstruct as a single microstrip track due to their 
extremely small opening angles, but will reconstruct as two separate MWPC 
tracks as the e"^ and e~ are bent in opposite directions in Ml. 
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Track segments from the microstrip tracks and the MWPC tracks which fail 
to hnk are used to form Vees, Kinks, H~'s, and 0~'s. 



3.4 Primary Vertex 

In the early stages of data reconstruction a basic vertex algorithm is imple- 
mented using only the microstrip track information. These vertices are used 
for reconstructing different Vee categories and are meant to be roughly cor- 
rect. Later in our analyses the "true" primary vertex is recalculated using all 
available tracking information with multiple scattering information included. 
The exact algorithm depends on the final state that is being investigated. For 
instance we use a candidate driven approach when the decay daughters are 
fully reconstructed such as D° — > K~7r~^, a special seed plane approach for one 
prong final states such as D'^ — > K'^n'^, and a stand alone vertexing algorithm 
for scmilcptonic decays such as K~7r~^ fi'^u. These techniques have been 

described in several E687 theses. A good reference is the thesis [3] of Luca 
Cinquini. 

The basic vertex algorithm minimizes the distance of closest approach of the 
tracks in the transverse plane. Specifically, the minimized quantity is: 

^2 ^ j^/ X-{x, + 4Z) y ^ ^ Y-{y, + b[Z) 

where X, Y, Z arc the coordinates of the vertex taken from the parameters of 
the fit, a^, b^, Xi, yi are the slopes and intercepts of the i-th track, and ax,i, o-y,i 
are the errors returned by the purely geometric track fit to the hits of the i-th 
track. 

The algorithm begins by assigning all microstrip tracks to a common vertex 
and by computing the corresponding x^- Tracks are removed from this vertex 
one at a time, beginning with the one that gives the largest contribution to the 
X^- Each time the vertex is refit and the subtraction process continues until 
the f^-lls below the threshold. Since the space location may have changed 
in the process, all discarded tracks are individually tested to check whether 
they originate from the found vertex and if so, they are included. Once the 
construction of the first vertex is completed, the procedure is repeated with 
the remaining set of tracks. At the end of the process each microstrip track is 
assigned to just one vertex, or they might remain unassigned. For Vee finding 
algorithms the most upstream microstrip vertex within the target is called the 
primary vertex. If no primary vertex is found, then the center of the target is 
selected as the primary vertex. 




9 



4 Vees 



Kg and A° (usually referred to as "Vees") are often found among the decay 
products of charmed particles. In FOCUS, these particles are reconstructed 
through the charged decay modes [4] : 

K'^s 7r+7r-(BR = 68.6%) 

AO 7r-(BR = 63.9%) 

These particles are relatively long-lived with respect to charmed particles, and 
may travel several meters in the spectrometer before decaying. Depending on 
the region of decay, they leave topologically distinct tracks, but they must be 
reconstructed with different algorithms. In all, Vees can be reconstructed in 
the FOCUS spectrometer over a decay length of about 10 meters. A sketch of 
the regions where different algorithms are employed for Vee reconstruction is 
presented in Fig. 3. 




Fig. 3. A schematic drawing of the regions of the spectrometer in which Vees are 

reconstructed by different algorithms. The solid lines are the two-prong daughter 
products of the Vees. The dashed lines indicate that the parent particle originated 
from the target region. Note that target box (TGT) is simply for illustration. In the 
experiment it is segmented and interleaved with silicon planes. 

All Vee-reconstruction codes have in common the search for a pair of oppositely 
charged tracks which originate from a common point in space, the Vee decay 
vertex. The invariant mass of the pair is computed by first assigning the pion 
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mass to both tracks to test the Kg hypothesis. Next, the proton mass is 
assigned to the faster particle and the pion mass to the slower particle to 
check for a A". Initially, there is no Cerenkov identification applied, and the 
Vee requirements are intentionally left loose to allow for different degrees of 
signal to noise in individual analyses. 

While not all Vee categories were used in FOCUS analyses, we did develop 
algorithms to identify and to reconstruct Vees in all regions described in this 
section. The SSD Vees, Ml Vees, and One-link SSD Vees were used in several 
analyses. The RECON Vee and P34 Vee categories were used in E687, but 
were not implemented in FOCUS. This resulted from a reduction in the beam 
energy such that long-lived decays became relatively less important. It is also 
the result of running the beam at much higher intensity which led to increased 
noise in the chambers which in turn made the rate of false RECON Vees high 
and the algorithm less efficient. For a similar reason the MIC Vee algorithm 
was implemented in E687, but not in FOCUS. All MIC Vees are first found 
as Ml Vees and it required excessive computer time to search over all the 
remaining microstrip hits. Finally, the Single-linked MPWC Vees were only 
used in searching for E" and D," decays. 

4.1 SSD Vees 

SSD Vees are reconstructed from pairs of oppositely charged linked SSD- 
MPWC tracks originating from a common vertex. As the decays in this cate- 
gory occur close to the primary vertex, they tend to be lower momentum and 
accordingly have an excellent mass resolution. These decays are principally 
Kg and A° which decay upstream of the second SSD station. The secondary 
Vee vertex is required to He downstream of the reconstructed SSD primary 
vertex. Basically, these Vees are found in the same way [5] that we find charm 
particles. As such if matching hits are found in the target silicon planes then 
they are added into the sihcon track definition. The SSD Vees are the cleanest 
and the best defined category of Vees: the Vee track for this category has a 
resolution comparable to that of two combined SSD tracks. The standard de- 
viations for the Kg and A° mass distributions in this category are 3.6 MeV/c^ 
and 1.6 MeV/c^ respectively. 

4.2 Ml Vees 

The Ml Vees are composed of Kg and A° which decay between the last SSD 
plane and the first MPWC station, PO. They are reconstructed with pairs of 
unlinked MPWC tracks and are divided in three subcategories, according to 
the nature of their components: track-track, track-stub, and stub-stub. 
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The reconstruction algorithm is substantially the same for the three topologies. 
For each candidate pair of unlinked MPWC tracks, the intersection in the xz 
plane (non-bend) is first found; an iterative procedure then traces the two 
prongs through the Ml field and determines the y location of the Vee vertex. 
In the case of a track-stub Vee, the tracing also allows the computation of 
the unknown momentum and charge of the stub prong. For a stub-stub Vee, 
it is necessary to further constrain the Vee vector to originate from the SSD 
primary vertex and then the unknown momenta can be computed. Finally, a 
global fit using the full covariance matrices of the tracks (including multiple 
Coulomb scattering effects) is applied to each Vee candidate to provide a 
better estimate of the Vee decay vertex and the Vee momentum. 

Ml Vees are the most copious Vee category, accounting for over 70% of the 
total reconstructed Vee sample. Their mass resolution and vertex resolution 
are not as good as for the SSD Vees and are a very strong function of the 
angle that the Ml field makes with respect to the normal of the Vee decay 
plane. The resultant mass distributions often have very non-Gaussian tails. 
To correct for this effect, Vees are retained based on a normalized mass cut or 
the difference between the reconstructed and nominal Vee mass divided by the 
anticipated resolution for a given topology. The normalized mass distribution 
is much closer to a true Gaussian distribution. Our selection criteria for the 
Ml Vees is that the absolute value of the normalized mass be less than 5. 
The mass resolution for the Ml Vee category varies from 6.6 MeV/c^ for 
stub-stub Vees to 5.9 MeV/c^ for track-track Yees. 



4.3 RECON Vees 



RECON Vees are Vees which decayed between PO and P2. Because their decay 
region is further downstream, they tend to be a Vee category with higher 
momentum. 

RECON Vees are reconstructed using hits in PI, P2, P3, and P4 which have 
not already been used by the general MPWC pattern recognition (which only 
finds tracks with hits in PO). This requirement significantly speeds up the 
algorithm, but it reduces the efficiency in higher multiplicity events. First, 
track projections in the xz (non-bend) view are constructed and checked two at 
a time for intersecting between PO and P2. Projections which do not intersect 
in the desired region with any other projection view are disregarded. Then, 
projections in the U-view, V-view, and Y-view are formed and matched to 
the xz plane projections to form tracks in space, with only loose requirements 
on the x^/DOF of the track. Several track topologies are allowed: P1234 {i.e. 
tracks with hits in PI, P2, P3, and P4), P123, P234, and P23. Finally, tracks 
are combined pairwise and a global fit to the Vee hypothesis is performed. The 
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parameters of the fit are the five parameters for each track {x and y slopes 
and intercepts and the y bend angle in M2) plus the coordinates of the Vee 
vertex. In the case of P23 candidates (for which prongs are defined by a single 
point on each side of the magnet) it is also necessary to assume that the Vee 
originates from the SSD primary vertex in the target. RECON Vees sharing 
the same xz projections are arbitrated on the basis of their x^/DOF. 

This type of Vee is not used in FOCUS analyses. 
44 P34 Vees 

P34 Vees are Vees which decayed between P2 and P3 in the magnetic field of 
M2. Because their decay region is the most downstream considered, they are 
the Vee category with the highest momentum. 

P34 Vees are reconstructed using hits in P3 and P4 which have not already 
been used both by the general MPWC pattern reconstruction and by the hits 
used in the RECON Vee category. Because the MPWC hits can be erroneously 
pre-assigned to another track category, the reconstruction efficiency for P34 
Vees is not high. Track projections are formed in the xz (non-bend) plane and 
checked two at a time for intersecting between P2 and P3. Projections which 
don't intersect in this region are discarded. The projections from the U-view, 
V-view, and Y-view are formed and matched to the X projections to form 
tracks. In order to determine the unknown momenta of the two station tracks, 
it is necessary to further constrain the Vee vector to originate from the SSD 
primary vertex. 

This type of Vee is not used in FOCUS analyses. 
4.5 One-link SSD Vees 

One-link SSD Vees are reconstructed from the combination of a finked SSD 
and a MPWC track with an unlinked SSD track. These are Kg and A° which 

decayed before the second SSD station, but for which one of the decay prongs 
falls out of the Ml geometrical acceptance and therefore is not found in the 
PWC. 

The finked SSD to MPWC track and the unhnked SSD track are required to 

make a good space vertex with a confidence level greater than 1%. A primary 
vertex is used as a constraint. The significance of separation between the 
primary and Vee decay vertex (L/ctl, is computed and it is required to have 
L/(7l> 10. Also the primary-secondary vector is required to lie in the plane 
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of the two decay prongs. By knowing the momentum of the hnked track, the 
primary-secondary direction and the unhnked track direction, it is possible to 
balance the transverse momentum and compute the total momentum of the 
unhnked track (there is no two-fold ambiguity in the kinematics) . Finally, the 
invariant mass of the two tracks is computed for the Kg hypothesis and for the 
A° hypothesis. In the Kg case no Cerenkov requirements in the two prongs are 
imposed, while in the A" case the linked prong is identified by the Cerenkov 
algorithm [6] to have a light pattern more consistent with a proton hypothesis 
than a pion hypothesis. The A° candidates where the pion is linked and the 
proton is unlinked are not reconstructed, since this kind of Vee contained very 
little signal over an overwhelming background. These Vees have a Kg mass 
resolution of 4.7 MeV/c^. 



4.6 MIC Vees 



MIC Vees have the decay vertex between the second and the fourth (last) SSD 
station. The reconstruction algorithm starts by projecting unlinked MPWC 
tracks backward onto the SSD detector and by searching for unused hits in 
the last two stations. If one or more matching triplets of hits are found, the 
parameters of the track are recomputed by a global fit which uses both the SSD 
and the MPWC hits. The new reconstructed tracks are then checked two at a 
time to see if they originate from the same vertex in space, and a cut on the 
distance of closest approach (DCA) is imposed. Candidate Vees sharing one 
prong are arbitrated on the basis of minimum DCA. The Kg mass resolution 
for this category is 4.4 MeV/c^. 

This type of Vee is not used in FOCUS analyses. 



4.7 Single-linked MPWC Vees 



Single-hnked MPWC Vees share the same decay region of the Ml Vees (be- 
tween the SSD detector and PO), but are composed of one linked SSD to 
MPWC track and one unlinked MPWC track. The reconstruction algorithm 
is essentially the same as for the Ml Vees, and they are divided into three 
subcategories: track-track, track-stub, and stub-stub. For the Kg reconstruc- 
tion there was too much background in these categories and they were never 
used. However, for the A'' they were very important particularly when the the 
AO was a decay product of a 5 or a O . 
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4-8 Arbitration of Vee types 



In high multiphcity events it is possible for a track to be found in more than 
one Vee category This occurs when one or both of the daughter tracks in a 
Vee arc found in other Vcc candidates. The Vee candidates are arbitrated so 
that a single MPWC track is used just once in forming a Vee. In order for a 
Vee to be arbitrated it must first have passed the basic requirements such as 
being within the mass windows and satisfying the normalized mass cuts and 
it must have a track shared with another Vee that satisfies the same basic 
conditions. The decision to keep a Vee is made with the following criteria. 

• Two SSD Vees sharing a microstrip segment are arbitrated based on which 
Vee has the smaller distance of closest approach at the decay vertex location. 

• Two Ml track-track Vees sharing a MPWC track arc arbitrated based on 
which Vee has the smaller distance of closest approach at the decay vertex 
location. 

• Ml track-track Vees are always selected over Ml track-stub Vees. 

• Two Ml track-stub Vees sharing a MPWC track are arbitrated based on 
which Vee has the smaller x^/DOF of the fit. 

• Ml track-stub Vees are always selected over Ml stub-stub Vees. 

• Two Ml stub-stub Vees sharing a MPWC track are arbitrated based on 
which Vee has the smaller x^/DOF of the fit. 

4.9 Vee Refit 

Vees are reconstructed with only a rough knowledge of their actual production 
vertex. For the purposes of track finding the assumption is made that the 
tracks originate at the primary vertex (see Fig. 4). When a secondary charmed 
particle vertex is identified, then the Vee is refit to the true production vertex. 
This refit shghtly improves the Vee production angle and the Vee mass. 

Six histograms of Kg Vee types are presented in Fig. 5. The histograms are 
all plotted on the same mass scale and were made from the same small data 
sample (approximately 0.5% of the total). The {Ml stub-stub) Vees have the 
most background and the poorest mass resolution. The {Ml track-track) and 
{Ml track-stub) Vees are the most common categories. The {MIC) Vees have 
the least background, but the yield is comparatively low. The {SSD) Vees 
have the best mass resolution, but can have considerable background since 
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RecpnstmctedDirei 



Primary Charm Decay Vertex Vee Decay Vertex 



Fig. 4. A sketch of the charmed particle decay K^tt~^ tt~^ tt~ where the 

TT+vr" is presented to display the need for a Vee refit. The is first recon- 
structed to originate from the primary vertex, but a refit is necessary when it is 
reassigned to come from a new vertex. 



the decay occurs close to the interaction region. The {One-link) Vees also 
have good mass resolution due to their very low momentum, but these Vees 
have considerable background. 





I 







0.5 

Ml Stub-Stub 



q.6 




0.5 

Ml Track-Track 



0.5 

Ml Track-Stub 




One-Link 



Fig. 5. Histograms of a sample of six categories of K^^s which arc used in the 
analysis. The dominant categories in charm decays are the {Ml track-track) and 
{Ml track-stub) Vees. The categories with the best resolution are MIC, SSD, and 
One-link Vees. 
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Nine histograms of A° Vee types are plotted in Fig. 6. The data represents 
about 0.5% of the total A*' yield. For the A*^ decays it is useful to retain the 
{Single-linked) categories, since S~'s and fl~^s decay to A°'s and occasionally 
the proton in the A° decay is erroneously linked with the S~ or Q~ track 
segment in the SSD. 






Fig. 6. Histograms of a sample of nine categories of A*^'s which are used in the 
FOCUS analyses. The final row of Single-linked A'''s are not used in charm analyses 
involving direct charm decays to A^'s, but account for almost 15% of the S~'s and 
r2~'s decays. 
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5 Kinks 



The term "Kink" refers to the topology where one charged particle decays 
to another charged particle and a neutral particle. The neutral particle is 
undetected in the tracking detectors. 

A schematic of a Kink decay is given in Fig. 7. The parent particle with mea- 
sured direction cosines, a, /3, 7 and with mass, nip, decays to particles mi and 
m2. Particle mi has its momentum and direction cosines measured. Parti- 
cle 2 is neutral and goes undetected in the spectrometer. Through kinematic 
constraints and by assuming its mass the parent momentum is calculated to 
within a twofold ambiguity. 

Daughter Energy E j 

Reconstructed Daughter Momentum = p j(app bpp cp^) 
Daughter Mass = rrij 



Parent Energy = E ^^^^ 
Parent Momentum = P (ccP,pP, yP) 

Parent Mass = mp j/^ Daughter Energy = E 2 

\ Missing Daughter Momentum = p 
\ Assumed Daughter Mass = ni 2 



Fig. 7. A sketch of the Kink decay process where a parent particle decays to two 
tracks where one track is completely measured and the other track information is 
unknown. The direction of the parent track is known, but its momentum is unknown. 
Through kinematics the momentum of the parent particle can be calculuated to 
within a twofold ambiguity. 

The momentum of p2 is determined through conservation of momentum, 

P2 = {aP - api, f3P - bpi, -fP - cpi) (2) 

and the square of the energy of E2 is given by 

Ei = {E-Eif=pl + ml (3) 

Substituting the expression for p| we obtain 

{ml + mj-ml) + 2Ppi {aa + pb + 7c) = 2EEi (4) 
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Squaring this equation, substituting = E"^ — nip, and solving for P, we find 

MC ± ^M2C2 - (1 - C2)(m2 - M2) 



where M = ^^^^I^J-^, C = ^^^°"+f+"^\ and P > 0. Thus, we find the twofold 
ambiguity. 

The three decays S+ pn^, S+ nvr^, and S~ nir" are reconstructed us- 
ing the Kink algorithm. It should be noted that the T,~^{suu) and the T,~{sdd) 
are not charge conjugate partners and that the 8 MeV/c^ difference in their 
masses is important to include in the algorithm. The E"*" decays weakly to p7r° 
51.6% of the time and to n7r~^ 48.4% of the time. The E~ decays to nvr" es- 
sentiaUy 100% of the time. A complete hst of the Kink candidates considered 
is given in Table 2. The fraction of the meson decays that we would recover 
from the Kink algorithm is insignificant compared to the prevalence of the 
non-decaying ones. We therefore concentrated on the hyperon decays, E's, S's, 
and fi's. Initially, we concentrated on the E decays. The decays E," — > Att", 
fl~ AK~, and fl~ 5°7r~ will be further considered in Section 6. 

Table 2 

Decay topologies considered in the Kink algorithm. 



Type 


Decay 




or (cm) 


Comments 


1 


K- 




371.3 


muon identified 


2 


E+ ^ 


n TT"*" 


2.396 


neutron energy found 










in hadron calorimeter 


3 


K- 


tt-ttO 


371.3 




4 






2.396 


p identified by Cerenkov 


5 


s- ^ 


n 7r~ 


4.434 


neutron energy found 










in hadron calorimeter 


6 




li~v 


780.45 


muon identified 


7 




Avr" 


4.91 




8 




KK- 


2.46 


K identified by Cerenkov 


9 






2.46 





The Kink algorithm begins by looping over each unlinked microstrip track 
that points into the Ml aperture and by pairing it with every unlinked MWPC 
track which also points into the Ml aperture. MWPC tracks used in Vee can- 
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didatcs arc not considered. A preliminary Kink vertex location is determined 
by intersecting the microstrip and MWPC tracks in the xz plane. This choice 
provides a rough estimate of the x and z location, which is required to be 
downstream of the last microstrip station and upstream of PO. 

If the MWPC track passes through all five MPWCs, then its momentum is 
already determined from M2 and it can be traced through the field of Ml to 
the estimated Kink z position. If the Kink z coordinate is upstream of Ml, 
then the y distance between the projected 5-chamber track and the microstrip 
track is compared and required to be less than 2.5 mm. This helps eliminate 
spurious Kink candidates. The parent momentum can only be calculated by 
making a particular decay hypothesis and by solving the kinematic equations. 
This involves assuming the parent (S^) and daughter masses (including the 
missing neutral daughter) and balancing the momentum transverse to the 
parent direction. This results in a two-fold ambiguity in the parent momentum. 

If the Kink z position is within Ml, then the ambiguity can be broken and 
a unique solution is found. The MWPC track is traced upstream to the z 
of the Kink. The x and y coordinates of the Kink are then given by the 
traced position of the track at this z. The microstrip track is then traced 
downstream through Ml to the Kink position. The trace is iterated several 
times, using trial momenta, until a momentum is found which best traces the 
microstrip track to the Kink vertex. The transverse momentum is balanced and 
the kinematic equations arc solved. If there are two physical solutions for the 
parent momentum, then the one nearest to the momentum calculated by the 
iterative trace is used. In theory, one could avoid solving the Kink kinematic 
equations for decays in the magnetic field by simply using the microstrip trace 
momentum, however often very little magnetic field is traced through and the 
resolution on the microstrip track curvature (momentum) is not well defined. 
For these reasons the kinematic solutions were always used to determine the 
momentum. 

For 3-chamber MWPC tracks we consider only the case where the Kink posi- 
tion is upstream of Ml. The x and y positions of the microstrip track at the 
previously estimated Kink z position are used as the Kink vertex location. 
An iterative trace procedure, similar to the one discussed above, is used to 
find the best momentum for the 3-chamber track. Unfortunately, a two-fold 
ambiguity exists in the parent momentum calculation. 

All putative Kinks are reconstructed according to all three decay hypotheses 
(types 2, 4, and 5) from Table 2. Some hypotheses are rejected on the basis of 
particle identification of the charged daughter or on the basis of calorimetry in- 
formation. Charged particle identification is made by re-running the Cerenkov 
code with the momentum determined from the Kink reconstruction in order to 
obtain a new Cerenkov identification estimate {i.e. the momentum changes). 
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For the S"*" pn^ hypothesis the Cerenkov hght pattern is required to be 
more consistent with a proton assumption than the pion assumption by a fac- 
tor of 7. For the S"*" nn~^ and S~ nvr" hypotheses the pion is required to 
not be consistent with being an electron, kaon, or proton and the neutron must 
impact the hadron calorimeter. The neutron must deposit sufficient energy, E, 
surrounding the point where the neutron is expected to strike the calorimeters 
and must satisfy 0.3 < E/p < 2.0, where p is the momentum of the neutron 
calculated from the Kink kinematic equations. A detailed description of the 
FOCUS hadron calorimeter and its performance is found in Reference 7. 

There is no method to plot the mass of the reconstructed Kinks directly be- 
cause of the assumed mass constraint in the algorithm. To show the success of 
our algorithm we present the invariant mass plots in Fig. 8 for S~7r"'"7r"'" 
where the S~ is reconstructed through the n7i~ channel and for A+ S+tt+tt^ 
where the S"*" is reconstructed through the pir^ channel and through the mr^ 
channel. Care has been taken to weight the cases where there are two solutions 
such that the yields in the caption for Fig. 8 are correct. 




2.2 2.4 2.6 2.2 2.4 2.6 2.2 2.4 2.6 



GeV/c^ GeV/c^ GeV/c^ 

Fig. 8. Three invariant mass plots for the decays (a) S+vr+vr^ where 

S"^ —>■ piT^ with a yield of 915±50 events, (b) — > S+vr^vr^ where mT~^ 
with a yield of 854zb65 events, and for (c) A^ — > S^vr+vr^ where S~ mr~ with 
a yield of 654±42 events. 



6 S 's and 's 

The — i- A^n^ and — > A^K^ decays are reconstructed via several tech- 
niques. The S~'s and fl~'s decaying upstream of the silicon microstrip de- 
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tector arc reconstructed differently than those decaying downstream of the 
microstrip. If a Vec is not found, or if the A''— track vertex is located down- 
stream of the Vee vertex (i.e. due to poor vertex resolution in the event), 
then we use a technique referred to as 'multivees' where there are three un- 
linked tracks. Finally, if no fully reconstructed S~ or multivee is found, then 
we use a Kink algorithm to find S~ A°7r~ where the A° is unidentified. 
This algorithm is needed to reconstruct the H~ when A° m:^. Each of these 
techniques will be described below. 

The fully reconstructed algorithm uses a common set of requirements to select 
Vees as A°'s. The Vee daughter track with the highest momentum is considered 
the proton for the A° hypothesis and the Cerenkov algorithm is executed using 
momenta of the Vee tracks as determined by the Vee algorithm. 

6.1 Upstream Reconstructed S""'s and Q~'s 

Upstream S~'s and fi~'s are those which decay upstream of the first microstrip 
station (i.e. within the target or between the target and the microstrip detec- 
tor). They are also referred to as 'Type 1' decays. A schematic of a typical 
decay in the category is presented in Fig. 9. 

Linked MWPC tracks are paired with each A° Vee which satisfies the A° 
hypothesis. Upstream decays are reconstructed by intersecting the A° vector 
and the MWPC track and by demanding that the confidence level of this vertex 
be greater than 1% and that the be consistent with originating from 

a production vertex further upstream. The distance between the production 
vertex and the E'~/Q~ decay vertex is defined as L. The track is assigned the 
pion mass to form the invariant mass for S~ hypothesis and is assigned the 
kaon mass for the Q~ hypothesis. There are two requirements which are used 
to significantly improve the signal-to- noise for these S^/fi^ signals. The first 
is the significance of separation of the E~ /fl^ decay vertex from its production 
vertex, and the second is the isolation of the E~ /ft" decay vertex from other 
tracks. 



PO PI P2 





SSE 




Ml 


























IT 

- IT 



p 



Fig. 9. A schematic drawing in the bend view of the spectrometer of a E~ decay 
which occurs upstream of the silicon strip detector (SSD). Only the front part of 
the spectrometer is displayed. 
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The significance of separation of the 'E~ decay vertex from its production 
vertex is defined as the quantity L/a^, where is the error on L. By increas- 
ing the cut on L/ul {i.e. L/a^ >1, 2, 3, etc.) one obtains a cleaner sample. 
The basic algorithm requires L/az > 0. 

The isolation of the decay vertex from other tracks is tested by at- 

tempting to place other tracks in the vertex and by refitting it. A cleaner 
signal is obtained by requiring that the confidence level from the fit to the 
new vertex be less than a certain value. Invariant mass plots for upstream de- 
cays for the Att" and AK~ combinations are presented in Fig. 10. Note that 
there is considerable background under the Q~ signal which is almost entirely 
due to pion particle misidentification from the S~ Att^ decays. One should 
also note that fl~ yield is about a factor of 20 less than the S~ yield. 




1.3 1.305 1.31 1.315 1.32 1.325 1.33 1.335 1.34 

GeV/c^ 

S" Upstream SSD (Typel) 



1.66 1.665 1.67 1.675 



1.685 1.69 1.695 

GeV/c^ 



Q' Upstream SSD (Typel) 



Fig. 10. The invariant mass plots for the Avr and the AK combinations for the 
category where the decays occurs in front of the SSD detector. The plots represent 
the full FOCUS data sample. 



6.2 Downstream Reconstructed 5 's and ft 's 



Downstream S~'s and Q~'s are those which decay downstream of the last 
microstrip plane and upstream of the first MWPC plane. They are also referred 
to as 'Type 2' decays. A schematic of a typical S~ decay in this category is 
presented in Fig. 11. The decay distance for this category is more than three 
meters along the beam direction. A very important advantage to these decays 
is that the leaves a track in the microstrip detector before decaying. 

This track can in turn be used for finding charmed particle baryon decay 
vertices. 

The reconstruction algorithm begins by pairing each A° with every unlinked 
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Fig. 11. A schematic drawing in the bend view of the spectrometer of a H~ decay 
which occurs downstream of the sihcon strip detector (SSD). Only the front part of 
the spectrometer is displayed. 

MWPC track in the event. An estimate is made for the A'^-track vertex by 
computing the z intersection of the A" vector and the track in the non-bend 
view {xz) plane. If this vertex is downstream of PO or more than 50 cm up- 
stream of the target, then the combination is rejected. An estimate for the 
X and y positions of the vertex is given by the x and y of the A° vector at 
the given z of the vertex (because the A*^ is neutral its path is not deflected 
by the field of Ml). The unlinked track is next traced through the magnetic 
field of Ml to the position of the putative vertex. If the track is 3-chamber, 
then an iterative procedure is used whereby the track is assigned a different 
momentum for each iteration until a good trace is made to the given vertex. 
A better determination of the vertex z position is made by computing the 
distance of closest approach between the A° vector and track. (This iteration 
is important because the track can pass through only one side of the magnet 
making fringe field corrections more significant.) 

Next, the microstrip track of the candidate charged S^/Q^ is found. The sum 
of the momentum vectors of the A° and the track are used to form a candidate 
S~/Q~. Unlinked microstrip tracks are used and an attempt is made to match 
each one with the candidate E~/Q~. If the A° under consideration is a single- 
link Vcc, then the Vcc-linkcd microstrip is used. Each candidate microstrip 
track is traced downstream to the A°-track vertex; if the vertex is within the 
field of Ml, then the magnetic trace is used otherwise the microstrip vector 
is simply extrapolated to the vertex. Because the microstrip track has much 
better position resolution than does the Vee, a better vertex position can now 
be determined. The new z position is defined as the z where the microstrip 
track and the MWPC track make their closest approach. If the z position is 
downstream of PO or upstream of the target, the microstrip track is rejected 
as a candidate. Also, the vertex is required to be upstream of the A° 

vertex. 

To remove spurious matches in this algorithm, a cut is made on the candidate 
vertex. The candidate vertex is calculated in two ways: the point where the 
A*^ and the MWPC track make their closest approach, and the point where 
the microstrip track and the MWPC track make their closest approach. A cut 
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is made on the transverse distance (in the x-y plane) between these two puta- 
tive vertices. The second quantity on which a cut is apphed is the difference 
between the x and y slopes of the microstrip track and slopes given by the 
sum of the A° and the MWPC track momentum vectors. The x and y slopes 
must agree to within 4 milliradians of the momentum vectors. 

The MWPC track is assigned a pion mass for the S~ hypothesis and is assigned 
a kaon mass for the Vt~ hypothesis. For the Vt~ candidates, the MWPC track 
is required to be identified by the Cerenkov counters as being consistent with 
a kaon hypothesis. 

Invariant mass plots for downstream decays for the Att" and A.K~ combina- 
tions are presented in Fig. 12. In this region one finds that the Vt~ yield is a 
factor of 30 lower than the S~ yield. The ratio of the 0," yield to the S~" yield 
is reduced downstream of the SSD detector due to the shorter Q~ lifetime 
relative to the H~. 




Fig. 12. The invariant mass plots for the Avr" and the h.K~ combinations for the 
category where the decays occurs downstream of the SSD detector. The plots are 
for the full FOCUS data sample. 



6.3 Multivees 



A 'multivee' is composed of three unlinked MWPC tracks and one unlinked 
SSD track. A schematic drawing of a multivee decay is presented in Fig. 13. 
While the category was designed to select three prong decays such as K~ — > 
7r~7r~7r"'", it has proven to be useful in recovering S~ and Vt~ decays where 
the A" decay vertex is close to the decay vertex. It also works well in 
reconstructing Vees from decays which open in the vertical plane of the 
magnet. 
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Fig. 13. A schematic bend view drawing of the decay of the K tt+tt tt using 
three unlinked MWPC tracks and one unhnked SSD track. 

Initially, three unlinked MWPC tracks are intersected in the xz view and an 
unlinked SSD track with the closest distance of approach at the xz vertex 
is selected as a match. The algorithm is separated into decays which occur 
upstream of Ml and within Ml. 

If the decay vertex is upstream of the Ml magnet, then the SSD track is 
extrapolated in z to the vertex and each of the stubs are traced to the SSD y 
vertex position. If there are five chamber tracks, then they arc swum upstream 
and intersected with the SSD track to find a better vertex position. 

If the decay occurs within the magnetic field, then there must be at least one 
five chamber track. The vertex of the three unlinked MWPC tracks is found in 
the xz plane and the y location of the vertex is determined by swimming the 
tracks to the z location. If more than one five chamber track exists, then the 
z location is determined using the combined x + y information. All remaining 
unlinked stubs are traced to the yz vertex and their momentum is calculated. 
Next, the sum of the three unlinked track momenta is found and assigned 
to the microstrip track. This unlinked SSD track is now traced downstream 
to the z vertex. Successful candidates must be near the vertex of the three 
unlinked tracks and must agree in slope to within 4 milliradians in both the x 
and y views to the combined momentum vectors of the three tracks. 

In Fig. 14 the invariant mass distributions for 7r+7r~'"7r~, p7r~7r~, and pK~7i~ 
combinations arc presented. While the charge kaon decay events are not used 
in our analysis packages, the decays and decays are used. 



6.4 Kinks 



The algorithm for reconstructing the S~ — >■ A°7r~ decays as presented in Fig. 
15 is similar to the E decays discussed in Section 5. Unfortunately, there is 
no good technique to reduce the sample with additional constraints such as 
energy in the calorimeter or with Cerenkov cuts. For this reason we chose 
to only reconstruct events where the S~ decay occurred within the magnetic 
field of Ml. There are three advantages to using this subsample. First, the 
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Fig. 14. Invariant mass distributions for vr^Tr+TT", p7r~7r~, and pK~ Tr~ combina- 
tions. Signals for — > 7r"'"7r+7r~, E~ — > A^tt" where A" — >■ p7r~, and fl~ A^K~ 
where A'^ — > pvr" are clearly evident. 

reconstructed pion or kaon is a five chamber track and its momentum is well 
defined by its passage through M2. Second, there is no two-fold ambiguity in 
this category as the S~ also bends in the magnetic field. Third, the decays 
are well-separated from any material and the background from large multiple 
scatters is significantly reduced. 




PO PI P2 

Fig. 15. A schematic bend view drawing for a H~ decay inside of Ml where the A*^ 
is not reconstructed. 



6.5 E+ Mass Plots 



Invariant mass plots from the full FOCUS data sample for the charmed par- 
ticle decay S'tt+tt"'" for each of the four categories discussed above are 
presented in Fig. 16. The mass plots were found with a significance of sepa- 
ration cut of L/aL> A between the secondary and primary vertices. S^'s are 
selected such that there is no overlap between categories. Upstream S~'s or 
Type 1 decays have no overlap with the other categories because they decay 
before the SSD detector. Downstream S~'s or Type 2 decays are fully recon- 
structed and are the cleanest of the decays. If a Type 2 decay occurs in an 
event, the S~ Kink and multivee algorithms are not run. The next cleanest 
category is the S~ multivees. If a S~ multivee is found, then the S~ Kink 
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algorithm is not used. 




Fig. 16. Invariant mass plots for the charmed particle decay S+ — >■ S~7r"'"7r"'" using 
typel (upstream S~) decays, using type 2 (downstream E~) decays, using 
Kink decays, and using multivee decays. 

From an inspection of Fig. 16, it is clear that the Type 2 category dominates 
the signal with 246±20 events, followed by the Kink category (59±9 events), 
the Type 1 category (29±6 events), and the multivee category (25±4 events). 



7 Summary and Conclusions 

We have briefly described the tracking algorithms of FOCUS and how the 
high resolution silicon microstrip system is integrated with multiwire propor- 
tional chambers. Further, we have described the various techniques developed 
to reconstruct and A° decays in a multiparticle spectrometer. We have de- 
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scribed the two-fold ambiguity that occurs from S+ and S~ decays to a single 
charged particle and an unobserved neutral decay. Through these 'Kink' kine- 
matics we are able to observe charm baryon decays. Finally, we have combined 
the techniques of Vees and Kinks and used these techniques in the reconstruc- 
tion of S~ and Q~ decays. We believe the techniques decribed in this paper 
will prove useful to future experiments and may serve to show what is possible 
in multipurpose spectrometers. 
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